Inhomogeneous systematic signals in cosmic shear observations 
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We calculate the systematic errors in the weak gravitational lensing power spectrum which would 
be caused by spatially varying calibration (i.e., multiplicative) errors, such as might arise from 
uncorrected seeing or extinction variations. The systematic error is fully described by the angular 
two-point correlation function of the systematic in the case of the 2D lensing that we consider here. 
We investigate three specific cases: Gaussian, "patchy" and exponential correlation functions. In 
order to keep systematic errors below statistical errors in future LSST-like surveys, the spatial 
variation of calibration should not exceed 3% rms. This conclusion is independently true for all 
forms of correlation function we consider. The relative size the E- and B-mode power spectrum 
errors does, however, depend upon the form of the correlation function, indicating that one cannot 
repair the E-mode power spectrum systematics by means of the B-mode measurements. 
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I. MOTIVATION 



The power spectrum of the weak gravitational lensing 
distortions of background galaxies is quite directly re- 
lated to the power spectrum of intervening matter P, Q • 
The weak lensing (WL) power spectrum depends upon 
the linear and non-linear rates of growth of structure 
since recombination, and upon the redshift-distance re- 
lation produced by the expansion. These dependences, 
plus the straightforward theoretical framework, make 
WL a very attractive tool for the constraint of the post- 
recombination Universe, e.g. dark energy. Current 5- 
10% measurements of the WL power spectrum have al- 
ready begun to place interesting constraints [j| , and very 
much larger-scale projects are planned to reduce the sta- 
tistical errors on the WL signal to 1 part in 10 3 or lower. 

To reap the benefits of these large surveys, systematic 
errors must be well below the small expected statistical 
errors. WL measurements are subtle and difficult com- 
pared to most astronomical data analyses. There are no 
"standard lenses" on the sky, so calibration of the WL 
shear data is a significant worry. The finite point spread 
function (PSF) width tends to circularize the appear- 
ance of background galaxies, squelching the WL shear 
signal. This must be corrected analytically, and any er- 
rors in this process, or inaccuracies in the estimate of the 
PSF size, will lead to calibration errors. Huterer et al. 
P| investigate the effect of overall mean shear calibra- 
tion errors on cosmological parameter estimation. It is 
also likely, however, that there will be spatially varying 
calibration errors that are larger than the error in the 
mean calibration. For example, as the PSF size varies 
during a ground-based survey, the resolution parameter 

R = 1 §■ of the galaxies will vary (here a g is the an- 
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gular size of a target galaxy). If we fail to track this 
variation properly, the inferred shear will be modulated 
by a factor (1 — SR(0)/R(9)) 5j. In this paper we calcu- 
late the effect of spatially varying calibration errors upon 
the measured power spectrum, and determine criteria on 
these systematic errors which will have to be met if they 
are to be made negligible in future surveys. 

Other spatially varying errors could arise from photo- 
metric errors or Galactic extinction which will modulate 
the effective depth of the survey. Both effects may lead to 
errors in the (photometric) redshift estimation and source 
galaxy distribution, see e.g. 0, which would in turn lead 
to local modulation of the observed shear. While a mod- 
ulation of redshift depth is not strictly equivalent to a 
multiplicative modulation of shear, our calculations will 
still permit an estimate of the level at which depth mod- 
ulations become significant. 

The lensing shear of sources at some redshift z s is 
a 2-component tensor function of the angular variable 
0. As reviewed briefly below, the shear field can be di- 
vided into "E" and "B" modes corresponding to curl- 
free and divergence-free deflections, with corresponding 
power spectra Pe{1) and Pb(1) for stationary isotropic 
fields. Gravitational deflections, being derived from the 
scalar potential, will produce only E-mode power in the 
weak limit. It is therefore Pe{1) that will be used for 
cosmological constraints, with the B mode serving as 
the "canary in the coal mine" to alert us to potential 
non-gravitational sources of systematic error. A spa- 
tially varying scalar calibration factor will alter the am- 
plitude of the E-mode power and convert some into B- 
mode power. In ijllll we quantify this effect in terms of 
the 2-point statistics of the calibration systematic. In 
flIVI we present solutions using several models for the 
systematic error, and show how the deleterious effects 
are in general determined just by the rms amplitude and 
characteristic angular scale of the calibration errors. <|V] 
summarizes the results and the requirements upon future 
surveys that can be derived from these results. 

Some related calculations exist in the literature. 
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Schneider et al. |7f investigate the B-mode signal that 
is created by inhomogeneous source distributions, using 
a formalism similar to that employed here. We note that 
the inhomogeneous-source effect can be avoided by con- 
sidering only cross-correlations between source bins that 
are disjoint in rcdshift. The calibration inhomogeneity 
that we analyze here will likely not be so easily avoided. 
Vale et al. [8| conduct a numerical test of calibration 
inhomogeneity by modulating the shear seen in a ray- 
tracing simulation, and then calculating the resultant 
power spectra. We will test our analytic results against 
their simulated data. 

A rough target for calibration systematics is that their 
effect on Pe(1) be smaller than the expected statisti- 
cal errors. For a single-screen lens analysis, the uncer- 
tainty in Pe{1) averaged over an interval AlnZ = 1 is 
l^ 1 f s y.y^ 2 Pe{1) in the sample-variance limit. So for ambi- 
tious surveys with / s k y ^ 0.5, the power spectrum statis- 
tical errors are ~ 1 part in 10 3 at I = 1000. At higher I, 
the uncertainties due to shape noise and inaccuracies in 
the non- linear clustering theory will become important. 
The tolerances may be tighter when one examines the 
impact of power-spectrum tomography rather than just 
a single power spectrum. So a good goal is to have the 
calibration- induced error APe(1) be < 10 -4 Pe(7). 

II. SHEAR FIELD DECOMPOSITION INTO E 
AND B MODES 

Decomposition of a spin-2 field, such as shear or the 
Stokes parameters, into curl-free and divergence- free part 
was suggested to be useful for weak gravitational lensing 
studies by Stebbins 0, and for cosmic microwave back- 
ground (CMB) polarization by Kamionkowski et al. Q(J 
and Zaldarriaga and Seljak [ll|. Crittenden et al. [l2| 
and Schneider et al. study its use in revealing non- 
gravitational signals in weak lensing surveys. We briefly 
review the decomposition of the shear field into indepen- 
dent E and B modes, following the notation of Schneider 
et al. 0- 

The gravitational lens equation in the one-screen ap- 
proximation relates the detected direction 6 of photons 
on the sky to the (unobservable) direction f3 of photons 
emitted by a source: A = 9 — /3, where A is the de- 
flection angle scaled by a factor depending on the angu- 
lar diameter distances in the observer-lens-source system 
[l3|. The gradient of the deflection field, being a ten- 
sor of rank two, is usually decomposed locally into the 
trace, symmetric traceless part, and antisymmetric part 
as follows Ajj = k 5ij + 7^ + lo tij , where the shear ten- 
sor jij is symmetric and e»j is the Levi-Civita symbol in 
two dimensions. We denote partial derivatives with re- 
spect to directions in the tangent plane on the sky in a 
standard fashion by a comma. Thus we may express the 
convergence k and the rotation u> as linear combinations 
of derivatives of the deflection angle: 2k — A^i + A2,2, 
2u> = A1.2 — A24. Also, the shear components (71,72), 



defined as 71 = 711 = -722, 72 = 712 = 721, may be 
written as 271 = A11 — A2 i 2, 272 = A12 + ^2,1- More- 
over, we can write the deflection field as a sum of curl 
free and divergence free parts A = A + + A x which 
can be expressed as the gradient of a scalar potential (j>+ 
and the curl of a pseudoscalar potential 4> x respectively 
We designate as "E-mode" the curl- free deflection 
A + , which resembles an electric field pattern, and can be 
due to the mass distribution, ft produces the tangential 
shear pattern j + Q. On the other hand, the divergence- 
free "B-mode" deflection A x resembles a magnetic field 
pattern. This mode reveals in measurements as a "ra- 
dial" shear j x (i.e., 7+ rotated by 45°) and it cannot 
be generated by lensing in the single-screen approxima- 
tion. The potentials <fi+ and 4> x are closely related to the 
convergence k and rotation uj via the Poisson equation, 
V 2 0+ = 2 k and V 2 (j> x — 2w. Since the single-screen 
approximation is thought to be valid in cosmological sit- 
uations |14| , gravitational lensing information is confined 
to the E mode while the B mode should be zero. Thus the 
presence of non-zero B mode would be due to breaking of 
the single-screen approximation or, more importantly, to 
a variety of processes not related directly to lensing, such 
as measurement calibration errors 

lialialia, clustering 
of source galaxies 0, or their intrinsic alignments [Isf . 
For a more thorough discussion of E/B-mode decompo- 
sition see Crittenden et al. [l2| . 

In order to quantify the contribution of systematic un- 
certainties to the E and B mode power spectra we in- 
troduce a pair of two-point correlation functions £+(#) 
and £1(0), following Schneider et al. 0. They are linear 
combinations of correlation functions of E and B compo- 
nents of the shear, defined for each pair of galaxies with 
respect to the preferred coordinate system in which their 
positions are Q\ = (0,0) and 82 = (0,0): 

aid) = <7i(0i)7i(0 2 )) ± (72^1)72(^2)) . (1) 

Moreover, the correlation functions can be expressed 
as follows in terms of E and B-mode power spectra, Pe(1) 
and P B (l), defined as (k{1)k{1')) = (2tt) 2 6 u (l + I') P E (l) 
and (uj(1)lu{1')) = (2n f 5 D (l + I') P B (l): 

1 f°° 

&(9) = ^J o dl1 (Pe(1)±Pb(1)) Jo,4(W). (2) 

We can invert those relations and obtain power spectra 
expressed in terms of the correlation functions [in what 
follows we use a convention that upper sign in the sum 
on the right hand side refers to E-mode power spectrum, 
lower to B-mode]: 

/>oo 

p e : b(i) = « dee [e + (6)j Q (W)±e(8)Mm] ■ (3) 

Jo 

We do not consider cross-power spectrum of E and B 
modes as it will vanish due to parity conservation 0. 
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III. EFFECT OF SYSTEM ATICS ON E/B 
POWER SPECTRA 

Ideally, we would like to measure the shear field j(0) 
directly. What we observe, however, is the coherent el- 
lipticity induced on an ensemble of galaxies, which (a) 
is defined by the distortion g = 7/(1 — k); (b) is im- 
parted on galaxies that are not intrinsically circular, and 
(c) are viewed through a finite point-spread function 
(PSF). The measured shear field d(6) will in practice 
be modulated or contaminated by various observational 
effects [l4|. Although techniques for shear extraction 
from galaxy i mag es have been extensively developed and 
tested 0, [la, [jjj , there remain imperfections which can 
be detrimental to precision cosmology. 

Throughout the paper we assume that the observed 
field is related to the true shear field by a position- 
dependent multiplicative scalar factor 1 + e(6) such as 
will result from a misestimation of the "resolution" [f| or 
"shear polarizability" [l^. The systematic field e(8) is a 
random field assumed to have zero mean and described to 
the lowest interesting order by the two-point correlation 
function. Thus we express the observed field in terms of 
the shear and the systematics fields as 

d(0) = (l + e(0)) 7 (0). (4) 

This relation is local in real space, so it is going to couple 
modes of the shear field in the Fourier space, i.e. have 
some non-local effect on the relevant power spectra. We 
assume that the shear field 7 due to massive structures in 
the Universe is uncorrelated with systematics field e(0), 
which is a Galactic or instrumental foreground. The ob- 
served two-point correlation function £, d (9) can in this 
case be written as 

Z d (6) = (d(0)d(0 + 0)) (5) 
= (l + (e(0)e(0 + 0)» (7(0)7(0 + 0)) (6) 
= (l+e(6))?(6). (7) 

We have introduced two-point correlation functions £ 7 (9) 
for the shear field and £ e (9) for the systematics field. For 
simplicity we will assume that the systematics field e(9) 
is homogeneous and isotropic. In practice the assump- 
tion of isoptropy is not restrictive, as the effects of an 
anisotropic systematic could be approximated to first or- 
der by considering the azimuthally averaged correlation 
function. 

Correlation functions for the distortion field, £+(#) 
and £l(#), may be expressed as products of the corre- 
lation functions for the shear and systematics £±(9) — 
(1+C(d))t±{8) which follow from eqs. © and Q. 
We can rewrite eq. Q in terms of the distortion in- 
stead of the shear and then account for systematic sig- 
nals e(8). We split the observed E and B mode power 
spectra Pjf B (l) into two contributions as follows 

Pe,b( 1 ) = P IbV) + &PkB(l), (8) 



where the term P], B (l) is E mode (B mode) power spec- 
trum of the shear and AP| B (l) represents the contribu- 
tions to the E mode (B mode) power due to systematic 
signals. We focus on these error terms in the remainder 
of the paper. Using eqs. 10 and (JSJ) they can be written 
as 

poo 

API B (l) = 7T / d69?(9) [C + (9)J (19) ± C(0)Mm] ■ 
Jo 

We assume that the shear correlation functions £± re- 
ceive contribution from E mode only, i.e. P B (l) = P K (l), 
P B {1) = 0, since B-mode cosmological contributions are 
expected to be a few orders of magnitude smaller on 
scales > 1' 0- The systematic errors IS.Pe.bQ-) to E 
and B mode power spectra can be written as integrals 
over the convergence power spectrum P K {1) with a win- 
dow function We.b(^ q)- 

poo 

APe,s(0= / dqqP K {q)W E .B(l,q), (10) 
Jo 

where the window function depends solely on the corre- 
lation function £ e of the systematic modulation, and is 
given by 

W e ,b = \ [ d99C(9) [J (l9)J Q (q9) ± J 4 (l9)Mq9)} . 
z Jo 

(11) 

In the limit of a systematic that is completely correlated 
across the entire observation, i.e. a constant calibra- 
tion error, we have ^ c (^) = S 2 , where S 2 is the vari- 
ance of the calibration error. In this limit we obtain 
W E (l,q) = S 2 (? - 1 (5d(? - q) and W B (l,q) = 0, where we 
have used an integral relation for the Bessel functions 
J °° d9 9J n (l9)J n {q9) =q- 1 S D {q-l) HJ. Thus the er- 
ror contributions to E/B power spectra are AP B (l) = 
£ 2 P K (Z) and AP B (l) = in this case, and there is no 
conversion of E power to B power, as expected. 

Numerical simulations of calibration inhomogeneity in 
[8| are presented in terms of the aperture mass statistics 
M ap (P) and M X (R) with compensated filter defined in 
pH |22| . We produce analytic predictions for inhomoge- 
neous calibration errors for comparison with the numer- 
ical results of Q using the same filter as they did. 

IV. MODELING OF SYSTEMATICS 

We consider several potentially useful models of the 
correlation function of the systematic signal £ £ (#) and we 
examine the dependence of E and B mode power spectra 
(|10fl on the characteristics of £ £ (9) . The correlation func- 
tions considered here are analytically tractable and able 
to describe a wide variety of random processes leading to 
systematic signals. Each correlation function considered 
here is assumed to describe a stationary, isotropic random 
field. We assume that the systematic field has a finite 
variance E 2 . Equation (|11|) shows that AP b .b(1) oc E 2 
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if E 2 is a prefactor to some otherwise fixed functional 
form for £ e . Moreover, a correlation function in 2-D has 
to be bounded from below by a global minimum value 
of the Bessel function Jq(x). This condition is met by 
our models because they are assumed to be non-negative 

m 

We also introduce for each correlation function a char- 
acteristic scale R1/2 where the correlation function drops 
to 50% of its zero-lag value E 2 . 



A. Gaussian family 



is safely taken when the argument of Iq or I4 is greater 
than 1 or 100, respectively. Thus we obtain from (|14|l for 
small scales the following 



W E (l,q) 
W B (l,q) 



1 - 

31 



(17) 



The asymptotic expression l|17[) is useful when comput- 
ing the small-scale systematics contribution to E and B- 
mode power spectra (|10|) . Due to the exponential term 
in the window function l|17f) it is effectively a Dirac delta 
function 5y)(1 — q). Thus we can write eq. (|10|) as 



As a first model let us consider correlation function 
having a Gaussian shape with characteristic scale 9q 



AP E (l) 
AP B (l) 



E 2 P K (Z) 



31/)-2;-2 



(18) 



(12) 



We have R\j 2 — 2 9q In 2 for the Gaussian. The Gaussian 
is chosen because they are usually easy to handle analyt- 
ically. In this case an i nteg ral over scale 9 in eq. i|ll|) 
can be done analytically [24| and we obtain the following 
window function 



if we consider small scales compared to R\/2- The asymp- 
totic behavior of E and B-mode systematics power spec- 
tra is seen in fig. ^ It is notable that, for sufficiently 
large the systematic E-mode contribution APe is sim- 
ply a factor E 2 of the convergence power spectrum P K (l). 
Moreover, APe follows the shape of P K (l), whereas the 
B-mode contribution AP B drops rapidly. 



W E ,B(hQ) 



Y?9} 



e -hjlii 2 W) 



'lo(9 2 lq)±I 4 (9 2 lq)] 



where Io(x) and Ii[x) are the modified Bessel functions 
of the first kind of zeroth and fourth order respectively 
[20|. Because of the exponential growth of Iq(x) and 
I±{x) with x it is useful to rearrage terms in (|13|) and 
rewrite this equation as 



W e ,bM 



E 2 9? 



o e -ke 2 (i-q) 2 



h(0 2 lq)±h{9 2 lq) 



(14) 

where we have introduced functions I n (x) = e~ x I n (x). 

The large-scale amplitude of We,b(J, q) can be derived 
by noting the asymptotic behavior of the modified Bessel 
function for small arguments: Iq(x) ~ 1 and Ia(x) ~ 
x 4 /384 if x <C 1. Thus for scales large compared to R1/2 

when 6 l < 1, we obtain W E , B {l,q) ~ 5E 2 #g e~l? o /2 . 
When we consider power spectra APe, B (1) in this regime 
we get the following expression 



lz 2 9 



-VP K 11 = — 



APeMI) 



In the above we used the fact 
q #0 e _1 / 2e ° 9 has a maximum at q = 



dqqP K (q) e 



1 



1 fl 2 r? 



(15) 
(16) 



that the function 
1/6q and can be re- 



garded narrow around its maximum. The error we make 
using this approximation is less than 10% for R%/2 = 1° 
and 40% for R 1/2 =V. 

For small scales where 8q I 3> 1 and 9 q ^> 1 , we may 
use another asymptotic formula for modified Bessel func- 
tions which leads to I n (x) ~ \/2ttx~ 1 / 2 [2(J]. This limit 



B. Patchy 

Let us consider a systematic signal e(0) which is per- 
fectly correlated within circular patches of diameter 9a 
on the sky. This kind of systematics could arise if the 
survey is a mosaic of (circular) telescope pointings, and 
each pointing has a constant calibration error that is sta- 
tistically independent of all other pointings. For example, 
the impact of time- variable atmospheric seeing on galaxy 
shape measurements could produce such a pattern. As- 
suming such a model, we can compute the correlation 
function £ e (#) (independent of the specific distribution 
(pdf) of e(9) amplitude) 



C (9) = E 2 



2 

1 

7T 



E(0 o -t 



where H is the Heaviside step function and R1/2 ~ 
0.404#o- This type of correlation and its effect on E mode 
signal degradation and B mode generation were studied 
by [8J using ray-tracing simulations (their "sharp mod- 
ulation" model). Although they assumed square areas 
of correlation (the shape of CCD detectors), our ana- 
lytic model should match this well if we perform angular 
averaging over the square pattern to get an isotropic cor- 
relation function. 

In this case a closed form for the window function is not 
attained, so we have to rely on numerical integration. We 
can deal analytically with the window function (|11H in the 
limit of small scales Wo ^> 1 and q9o ^> 1. For this pur- 
pose we can approximate (fT§|) by £ £ (8) w E 2 [1 — 9/9o] 
where #0 = 2Ri/2- Let us use asymptotic formulae for 
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the Bessel functions for large arguments |20j and write 
the E mode window function as follows: 



W E (l,q) 

de 



1 — — | cos (w — ^- ) cos 
OqJ V 4 



-J- (20) 



Using the formulae for addition of cosines and subse- 
quently perform elementary integration leads us to 



W E (l,q) 
^1 - 



TT\/Tq 
cos 6q (I — q) 



91(1 -qf 



l-sin6> (l + q) 
2 S + q? 



(21) 



Thus in the interesting case of small scales main contri- 
bution to the window function comes from I w q which 
leads to 



W E (l,q) 



]20 Q sin 2 Wpl 



2nl ej(i- q y 



£ 2 r%(z-g). (22) 



The B-mode window tends to zero because of the identi- 
cal asymptotic behavior of Jo and J4. Thus in the small 
scales limit we have AP E (l) ~ E 2 P K (0, AP B (l) -> 
which was the case for a Gaussian correlations as well. 



C. Generalized exponential family 




FIG. 1: Statistical signal due to multiplicative shear errors. 
Shown are systematic-error contributions to the power spectra 
in E mode, APe(1), (solid) and B mode, APb(1), (dashed). 
The correlation function £ e of the multiplicative calibration 
field is assumed to be gaussian with 4% rms. The three plot- 
ted power spectra assume different characteristic scales R1/2' 
1°, 10' and 1'. At large scales, E mode and B mode contribu- 
tions are equal to each other for a given characteristic scale. 
Also shown is the convergence power spectrum P K (l) (dotted). 



A broad class of correlation functions can be described 
by a generalized exponential family [23j as follows 



e(0) 



l T(v) \8 



(23) 



where K v (x) is the modified Bessel function of the sec- 
ond kind, #0 is a characteristic scale and < v < 1. For 
v = 1/2 we obtain an exponential correlation function 
C(0) = Y?e- e ' e ° with R 1/2 = 9 In 2. When v < 1/2 the 
correlation function depends on 6 sub-exponentially on 
small scales and super-exponentially on large scales. For 
v > 1/2 the above behavior is reversed. Exponential- type 
correlations decay more slowly than do Gaussians with 
the same R1/2, so offer a test of the generality of the be- 
havior of AP e .b(1) for a given characteristic scale Ri/2- 
The window functions W EiE also decay slowly compared 
to the Gaussian case ( HIV Af) . 

For a generalized exponential family (|23|l we can com- 
pute analytically the respective power spectra 



P*(l) 



47rS 2 f 



1 + (io y 



-("+1) 



(24) 



The power spectrum has power-law scaling at small 
scales: P(l) oc I" with n = -2(v + 1) for W Q > 1. 
The allowed range of spectral indices is —4 < n < —2 



(a power-law correlation function must have n > —2). 
An example of a systematic signal of this type could be 
dust extinction in our Galaxy. Schlegel et al. [§5] show 
that the extinction pattern on the sky can roughly be de- 
scribed by the power-law power spectrum P(l) oc l~ 5 ^ 2 
for scales larger than ~ 15', corresponding to v = 1/4. 



V. RESULTS 

The cosmological background model we assume is 

ACDM with O ro = 0.3, fl A = 1 - Q m , H = 

TOkms^Mpc^ 1 , er 8 = 0.93. The distribution of 

source galaxies in redshift is assumed to be dN/dz oc 



z 2 exp 



(z/z 



,3/2 



with zq — 2/3 and mean redshift 

1. These are as assumed by Vale et al. Q, so that we 
may test our results against their ray-tracing simulation 
results. We compute the convergence power spectrum 
P K (l) using fitting formula for 3-D dark matter power 
spectrum given by Smith et al. |2rj| . 

In fig. 2]we show the power spectrum contributions due 
to an inhomogeneous calibration field described by the 
Gaussian correlation function i|12|) for three characteristic 
scales R\/2'- 1°, 10', and 1'. We take the rms of the 
systematics field to be S = 4%. Recall that AP EjE (l) oc 
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E 2 . We notice that APe spectra are featureless and have 
maxima near the maximum P K (l) (except for very small 
Ri/2)- On the other hand, the B-mode power spectra 
APb have maxima near the characteristic scale of the 
correlation function. 

In order to assess whether the signal due to systematics 
can be potentially harmful for weak lensing results, let 
us compare the contaminating power spectra APe,b(1) 
to the statistical errors on the convergence power spec- 
trum 5P K (l) [22]]. Assuming gaussianity of the conver- 
gence field we have with sufficient accuracy for our pur- 
pose that 



SP K (l) 



1 



(/ AZ/ sky ) 



1/2 



P K (l) 1 + 



n g P K (l) 



(25) 



where / s ky is the fraction of the sky covered by a survey, 
n g is the density of source galaxies with measured shapes, 
and ct 7 ~ 0.3 is galaxy shape noise. The P K (l) data will 
have to be binned over some interval Al for a mean- 
ingful comparison with the systematic error APe,b(J)- 
Because P K (l) is virtually featureless and there is no cos- 
mological information in its detailed structure, we choose 
broad bins of width Al — I. An even broader binning 
scheme would lower the 5P K (l) line in Figure [21 and our 
derived requirements on E would scale as (Al) -1 / 2 . Fu- 
ture, ground based, wide-field surveys like LSST [3l| are 
expected to cover / s k y ~ 50% of the sky and obtain good 
shape measurements for about 30 galaxies per arcmin 2 . 
Figure [5] shows the convergence power spectrum P K (l) 
and its statistical errors l|25|l for these values of / s k y and 
n g . The encouraging implication of fig. [21 is that keeping 
systematics (e.g, shear calibration errors) below 3% rms 
(E < 3%) should avoid significant contamination of the 
observed P K (l), even for future surveys. 

Figure [21 plots systematic-error power spectra 
APe,b(1) for a gaussian i|12|) . "patchy" l|19|l . and expo- 
nential (|23|l correlation functions, each with R1/2 = 10' 
and E = 4%. We notice that the shape of APe,b(1) is 
nearly independent of the specific shape of the correla- 
tion function e(6) of the systematic field. Thus from a 
practical point of view the important features of the sys- 
tematic field are the characteristic scale of correlations 
and the rms of the field. The latter affects the overall 
amplitude of the systematic errors as APe,b(1) oc E 2 . 
The former fixes the amplitude of the APe(1) — APb{1) 
at large scales, and gives the scale where the B mode 
starts decaying. 

We can compare our analytic results for the "patchy" 
correlation function to the numerical tests of Vale et al. 
[3. We set 6> = 25' and E = 10% to match the 
calibration-error pattern they superpose on their ray- 
tracing data. Our analytic estimates of the errors in- 
duced in the aperture mass variances (M 2 p x )(i?) are 
shown in fig. [21 These errors can be directly compared 
to those shown in Figure 2 of Vale et al. Q, which we 
reproduce in our Figure. Our estimates closely reproduce 
the results of the numerical simulation, except that we do 




FIG. 2: Statistical and systematic errors in the conver- 
gence power spectrum P K (l) (dotted): the statistical uncer- 
tainties SP K (l) are due to sample variance and source galaxy 
shape noise (thick solid). The systematic uncertainties in E 
mode, APe(1), (thin solid) and B mode, APs(i), (dashed) 
are shown for different calibration correlation models: Gaus- 
sian, "patchy" and exponential. They all share a characteris- 
tic scale R1/2 = 10' and E = 4% rms amplitude. Note that 
the E-mode systematic errors are essentially independent of 
the functional form, and are at worst equal to the statistical 
errors. 



not produce trough of A (M 2 p ) at the characteristic scale 
around 25 arcmin. The trough might be attributable to 
sample variance from the finite number (64) of patches 
used in the simulated images. 



VI. CONCLUSIONS 

We consider the effect of spatially varying multi- 
plicative systematic errors (assumed uncorrelated with 
the cosmological signal) on the measured power spectra 
Pe(1) and Pb(1) in the case of the 2D lensing. The prime 
example of this type of systematic would be shear cal- 
ibration errors which vary across the survey area due 
to changin g ob serving conditions. As shown by Hirata 
and Seljak |15| overall shear calibration errors of existing 
methods of shear measurement can reach 10% for galaxies 
of size comparable to the PSF. Such errors grow larger if 
one uses more poorly-resolved galaxies, as would be the 
case for deep ground-based surveys like LSST. Uncor- 
rected Galactic extinction could also introduce spatially 
correlated systematics in survey depth, altering the ob- 
served shear correlation functions. When we examine a 
variety of functional forms for the correlation function £ e 
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FIG. 3: Aperture mass variances as computed analytically 
in this work (solid) and obtained from ray tracing simula- 
tions [8j (dashed): there is a good agreement between these 
two approaches. The two upper curves, denoted by (Mf p ) 
represent the signal from the cosmological E mode (without 
systematics) and the lowest curve, (My), shows the cosmo- 
logical contribution to B mode from simulations (in our case 
this contribution is zero). The two pairs of curves denoted 
A (Mfp) and A(Mx) show the systematic signal contribu- 
tion to the aperture mass variance in E mode and B mode, 
respectively. The underlying modulation is of "patchy" type 
with 0o = 25' and 10% rms (sec. ITTBl . 



Comparison of the systematics errors APe(Z) on the 
power spectrum to the statistical errors expected for fu- 
ture weak lensing surveys indicates that we should not be 
afraid of systematic contamination if we keep calibration 
errors below E 3%. 

The absence of B-mode contamination in the most 
recent cosmic-shear measurements ^3 H3 suggests 
that systematic errors of the type considered here are 
below current statistical errors (5-10%) and hence do 
not bias the conclusions. Note, however, that B-mode 
power APb(1) consistent with zero on scales I ^> 1/0q 
does not necessarily imply the absence of significant cal- 
ibration error APe(1) (see Figure [TJ. Hence the present 
cosmic shear results could be significantly affected by cal- 
ibration errors, if they have a correlation length 9q that is 
larger than scales considered in the B-mode measurement 
(9q ^> 1//). Future surveys will beat down statistical er- 
rors, so we will have to understand and beat down sys- 
tematic errors as well. This work suggests that spatially- 
varying calibration errors will have to be reduced to 3%. 
This is well below the levels that have been demonstrated 
to date, but is probably achievable for well-behaved data 
with careful shape- measurement techniques |29ll30| . 
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